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The Chern-Simons approach has been widely used to explain fractional quantum Hall states in the framework 
of trial wave functions. In the present paper, we generalise the concept of Chern-Simons transformations to 
systems with any number of components (spin or pseudospin degrees of freedom), extending earlier results for 
systems with one or two components. We treat the density fluctuations by adding auxiliary gauge fields and 
appropriate constraints. The Hamiltonian is quadratic in these fields and hence can be treated as a harmonic 
oscillator Hamiltonian, with a ground state that is connected to the Halperin wave functions through the plasma 
analogy. We investigate conditions on the coefficients of the Chern-Simons transformation and on the filling 
factors under which our model is valid. Furthermore, we discuss several singular cases, associated with states 
with ferromagnetic properties. 

PACS numbers: 73.43.Cd, 71.10.Pm, 11.15.-q 



I. INTRODUCTION 

In understanding the fractional quantum Hall effect 
(FQHE), the now famous trial wave function proposed by 
Laughlin 1 proved to be a successful approach to describe the 
physics of incompressible quantum liquids at certain frac- 
tional filling factors. Laughlin's wave function is further- 
more the inevitable starting point for several generalisations, 
such as Jain's composite-fermion proposal, 2 '- Halperin's two- 
component wave function, 4 or more complicated wave func- 
tions describing states possessing quasi-particle excitations 
with nonabelian statistics^ 

A field-theoretical approach, complementary to the above- 
mentioned one, consists of so-called Chern-Simons theories, 
which formalise the idea of flux attachment that is also im- 
plicit in the trial wave functions. Chern-Simons theories have 
been successfully elaborated to study incompressible^ and 
compressible 7 quantum liquids in one-component systems as 
well as two-component systems, 8-12 which comprise e.g. bi- 
layer quantum Hall systems or single layer systems in situ- 
ations where the spins are not completely polarised. Multi- 
component Chern-Simons approaches have also been pro- 
posed in the study of edge excitations of the incompressible 
quantum Hall liquids. 13 An undeniable advantage of these 
Chern-Simons theories consists of their transparent insight 
into the exotic properties of these quantum liquids, such as 
their topological degeneracy, the fractional charges of their 
quasi-particle excitations or the statistical properties of the 
latter. 8,9 However, the Chern-Simons theories are usually less 
adapted when it comes to calculating quantities involving en- 
ergy scales. Indeed, Chern-Simons transformations act on the 
kinetic part of the electronic Hamiltonian, whereas they leave 
the interaction part invariant. The kinetic part gets therefore 
renormalised but continues to determine the overall energy 
scale whereas the physical energy scale in the FQHE must 
be set by the electron-electron interactions. 

A successful generalisation of Chern-Simons theories, that 
does not suffer from the problem of the correct energy scale, is 
the Hamiltonian theory proposed by Shankar and Murthy^— 
This theory is a very powerful tool for the computation of 



physical quantitiesr 7 -^ and even for the description of higher- 
generation composite fermion states. 19 However, it is limited 
by the fact that it does not incorporate internal degrees of free- 
dom. The success of the single-component Hamiltonian the- 
ory justifies a generalisation that can be applied to describe 
systems for which internal degrees of freedom (spin and/or 
pseudospin) are relevant. The main interest in such a gen- 
eralisation stems from realistic systems with more than two 
internal degrees of freedom, such as graphene with its four- 
fold spin-valley degeneracy 2 ^ or bilayer quantum Hall sys- 
tems with non-polarised electron spins. 

In this paper, we analyse a multi-component Chern-Simons 
theory within the framework of the microscopic theory by 
Shankar and Murthy. 14-18 This approach has two main ad- 
vantages over the previously proposed ones. First, it allows 
one to distinguish between physically relevant Chern-Simons 
theories from those which are ill-defined. The basic ingredi- 
ent for this distinction is the n x n charge matrix K, which 
was first introduced by Wen and Zee. 8,9 We find that ma- 
trices with negative eigenvalues need to be discarded in the 
study of physically relevant Chern-Simons theories because 
they would lead to ground-state wave functions that cannot 
be normalised. This structural feature of Chern-Simons theo- 
ries finds its physical interpretation within Laughlin's plasma 
analogy 1 that indicates a tendency of the different components 
to undergo a phase separation and thus to form spatially in- 
homogeneous states. We show that zero eigenvalues of the 
charge matrix K, in contrast to the unphysical negative eigen- 
values, find a compelling interpretation in terms of ferromag- 
netic quantum Hall states. Our results thus generalise previous 
work on two-component systems by Lopez and FradkiniS to 
an arbitrary number of components k. 

A second advantage of the present approach consists of 
a transparent connection between multi-component Chern- 
Simons theories with trial wave functions. It has been shown, 
in the simpler one-component case, that treating the fluctu- 
ations of the Chern-Simons vector potential within the har- 
monic approximation (Gaussian model) yields Laughlin's and 
Jain's (unprojected) composite-fermion wave functionsji^^i 
Similarly, we obtain here, within the Gaussian model of 



2 



K-component fluctuating Chern-Simons vector potentials, 
multi-component trial wave functions 22 that are generalisa- 
tions of Halperin's two-component wave functions4^2 Fur- 
thermore, we obtain in the same manner composite-fermion- 
type wave functions that may be viewed as particular multi- 
component generalisations of Jain's original proposal.^ 

The paper is organised as follows. In Sec. UU we define 
the Chern-Simons transformations for systems with k com- 
ponents and introduce, in Sec. [HI] extra degrees of freedom, in 
the form of the auxiliary gauge fields, as described by Shankar 
and Murthy. We subsequently diagonalise the harmonic os- 
cillator Hamiltonian and investigate the connection of the re- 
sulting wave function with trial wave functions through the 
plasma analogy. In Sec. lIVl we extend our results to the situ- 
ation of singular K matrices and discuss the relation between 
residual symmetries and underlying ferromagnetic properties 
of the quantum Hall states. Our conclusions are presented in 
Sec.|V] 



II. CHERN-SIMONS TRANSFORMATIONS 

We consider a quantum Hall system with k internal states, 
hereafter referred to as "components". In the simplest case of 
a two-dimensional electron gas at a GaAs/AlGaAs interface, 
one has k = 2 for the two possible orientations of the electron 
spin. The case k = 4 is relevant for bilayer quantum Hall sys- 
tems, where a second pseudospin mimics the layer index, or 
in graphene due to its two-fold valley degeneracy, in addition 
to the physical spin of the electrons. Higher values of k are 
rarely discussed in the literature, but may play a role in the 
context of multilayer systems or of bilayer graphene, where 
the zero-energy level consists of the n = and n = 1 Landau 
levels. 20 The Chern-Simons transformation 2 - 23 is defined by 
the relation between the n original electronic fields ip a (r) and 
the k transformed fields tp2 s (r) as 

^(r) = expf-i ( d 2 r'0(r-r')^^ Q ^(r'))^ S (r), 

(2.1) 

where 6(r) — arg(x + iy) indicates the angle between 
the vector r = (x,y) and the direction, and pp(r) = 
ipl(r)i/jp(r) = i/ip ^ (r)ijjp S (r) is the density operator of the 
particles of component /3. The n x n matrix K a p encodes 
the topological properties of the underlying quantum liquids, 
such as its degeneracy, the charges of its quasi-particle ex- 
citations and the statistics of the latter.^ Physically, it indi- 
cates the number of flux quanta attached to particles of com- 
ponent a due to the density of particles of component f3. This 
transformation is a singular transformation for the reason that 
6{r — r') has a singularity at r' = r. 

The gauge transformation is defined such that it generates 
the gauge potentials 2 ^ 




and such that the one-particle Hamiltonian [— ifiV + 
e A(r)] 2 /2m for the component a is transformed to 

H a = ±- 1-iKV + e A(r) + e A™(v)] 2 . 

Here, m is the mass of the particles, and e is the electron 
charge. By using V x V0(r) = 2ir S(r), we derive the corre- 
sponding magnetic fields, 

B° s (r)=--£tfa/ J ^(r)e a . 

6 P 

Since A^ s is a gauge field, its Fourier transform A^ s (q) 
may be fixed to a convenient gauge. We choose it to be trans- 
verse, iq-A^ s (q) = 0, so that it fixes the direction of A^ s (q) 
to be e z x q/|q|, up to a sign. For the magnitude, we use that 
under a Fourier transform B(r) = V x A(r) transforms to 
B(q) = iq x A(q), such that we obtain 

e iq| p 

where we define the transverse unit vector as e^r = ie z x 
q/|q|- 

The effective magnetic field seen by the composite particles 
of type a is 

B* a = B + (B^ s )=B{l-Y / K a 0^)e z , (2.4) 



where up are the component filling factors, given by v$ = 
(h/eB)ri0 — 2Trl B np in terms of the electronic densities 
ri0 and of the magnetic length Is = yKfeB. This re- 
sult is an extension of the two-component case presented in 
Refs. [TllfT2l Notice that each particle type has its own effec- 
tive magnetic field, and hence also its own magnetic length 
Ib* = y/h/eBa- The composite particle filling factors v* a are 
expressed in terms of the electronic filling factors v a asil 

12L = JLsl = _ = . . (2.5) 

v a 1% B* a 1-Y I K <X0V0 

This result generalises the one-component relation 

v* = <-> v = , (2.6) 

l-2sv 2sv* + l 

in terms of the Chern-Simons charge K = 2s. 

The statistical angle associated with the exchange of the 
transformed fields ip® s an ^ V'S 8 ^ can be derived by using 
their definition, Eq. d2.U . and the fact that the original fields 
are fermionic. Under the condition that the charge matrix K a p 
is symmetric, which is a generalisation of the condition dis- 
cussed in the two-component case, 11 we obtain 

^ S (ri)^ S (r 2 ) + e^^^f (r 2 )^ S (n) = 

and 

€ S (r i )^ S f (r 2 ) + K -' ^ S * (r 2 ) i/£ s (n ) 

= S a p S(ri - r 2 ). 
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Thus, we have found that the statistical angles of the exchange 
are nK a p, i.e., proportional to the entries of the charge ma- 
trix. The parity of the diagonal elements K aa of the charge 
matrix K determines the statistical properties of the Chern- 
Simons fields ip2 S - If mev are even integers, the originally 
fermionic electron fields tp a are transformed into fermionic 
Chern-Simons fields. However, one may also change the sta- 
tistical properties of the fields from fermions to bosons by us- 
ing odd integers for the diagonal components K aa . In the fol- 
lowing sections, we mainly discuss fermionic Chern-Simons 
fields, in order to make a connection with the composite- 
fermion theory, although the main conclusions of the paper 
also apply to bosonic fields. 



III. GAUSSIAN THEORY 
A. Auxiliary gauge fields 

The formalism proposed by Shankar and Murthy 1518 allows 
us to treat the fluctuations of the Chern-Simons vector po- 
tential via the introduction of n real-valued transverse gauge 
fields a° (r Jci The extended Chern-Simons Hamiltonian in 
first quantisation, with N a particles of each type a, reads 



H CS = 



a j a =l 



+eSA^(r ]a ) + e a ° a (r Ja )Y, (3.1) 

where we absorb the average value of the Chern-Simons 
potential ( I2.2l i into an effective vector potential A* (r) = 
A(r) + (A^ s ). This definition yields the effective magnetic 
field V x A* (r Ja ) = B* (r ,- a ) given in Eq. d2.41 >. In Fourier 
space, the fluctuations <5A^ S (q) are transverse, similar to the 
gauge field itself, as given by Eq. ( 12. 31 . Here, we have 



«SA£ S (q)=<$4f(q)e. 



cs. 



h 



K a p Spp{q)e A c 



Since we have artificially added the auxiliary gauge field 
a° (r), we have enlarged the Hilbert space, where the physical 
states form only a subspace {|</> p hys)} characterised by 



a°(q)l 



9 phys/ 



0. 



(3.2) 



for all components a. In other words, the gauge field operator 
acting on any physical state vanishes. 

Additionally, we introduce a longitudinal field P°(q) = 

iP°(q)e q (with eq = q/|q|), conjugate and perpendicular to 
the newly introduced gauge field a°(q) = a°(q)eq, accord- 
ing to the commutation relation in Fourier space 

K(q),P|(-q')] =»^V- 

Since the operator P° is conjugate to a° , it generates transla- 
tions in a° , as may be seen from the definition 



c/ = cx p(^EE^°(- c i') M S s (q') 

^ a q' 



which translates by the vector — 5A^ s (q) as U^a^ (q)U = 
a^(q) — 5A^ s (q). By using this shifting property of 
U, which is also valid in r-space, and with [p,- Q , U] = 
(h/eL 2 ) J2fs ^a^P^jJ, we may eliminate £A^ s (q) from 
the Hamiltonian of Eq. (13. IK which then transforms into 



1 N a 



Pj« +eA*(r,J 



<^J + ^E^ P ^) 



while transforming the states to ip CP = U^ 1 ^ 5 . In these 
equations, L 2 is the area of the system. By transforming the 
states, we also transform the constraint ( 13.2b to 

(a°(q)-^ s (q))|^ phys ) 

"„<ql - TTTtE^"/ 3 ^^ ) l^phys) = °- < 3 - 3 ) 



The Hamiltonian may be decomposed into three terms, 

H C p = H* + H coupi + P aux , given by 



H * =—yyu 



A 2 

2m ^ ^ 3 " 

a j a = l 



(3.4) 



H, 



coupl 



^EEX' ea^J + ^i^P^-J 

(3.5) 

i Na r 

^EE <W*J (3-6) 



b2 E E p l ^ ) K ^K a ,p° ( Tja ) 

(9 7 



where Uj a = pj a + eA* (r j a ) and b = h/ eL 2 , which has the 
dimensions of a magnetic field. Notice that for P aux we have 
used that a° (q) and (q) are perpendicular. 

In the remainder of this paper, we discuss only the term 
P aux that involves the auxiliary gauge fields. The full the- 
ory, including the other terms of the Hamiltonian shall be dis- 
cussed in a future publication. 25 



B. Gaussian model of the auxiliary gauge fields 

We shall now analyse P aux in detail. By observing that 
J2j a ^( r — r j Q ) = Pa{ r ) — n a + Sp a (r), we can rewrite 
Eq. (13.61 1 as 



+ ^^p;(r)^v;(r) 

1 J 
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Up to this point, all equations are exact. Now, we approxi- 
mate P aux by assuming that the density fluctuations Sp a are 
small with respect to the average densities n a . Since the re- 
sulting Hamiltonian becomes quadratic, this approximation is 
called the harmonic approximation. We should keep in mind 
that this approximation breaks down if the fluctuations are not 
small with respect to the average densities. In particular, the 
approximation is certainly invalid if one of the average densi- 
ties is zero. We therefore assume that none of the average den- 
sities n a vanishes. However, in the case of a singular charge 
matrix K, a redefinition of the filling factors might lift this 
problem, as will be discussed in more detail in Sec. [TV] The 
Hamiltonian P aux in Fourier space is approximated by 

tfo SC = ££^(e 2 <(-qK(q) 

q a ^ 

+ 52 EE^°(- c i)^^^ 7 (q)) I (3-7) 

where we note that o°(— q) = (a°(q))t and P°(-q) = 
(P°(q))T. Because the Hamiltonian (13.7b is quadratic in the 
gauge fields a° and its conjugate fields P°, it is possible to 
write it in terms of ladder operators. However, due to the 
appearance of the matrices K in the term with P°'s, it is 
a nontrivial task to define suitable ladder operators *4 Q (q) 
such that the commutators between them are of the form 

[A*(q),^3(q')] = $ a p$q,q>- 

In order to diagonalise the Hamiltonian, we define N = 
diag(jV Q }) as the dimensionless diagonal matrix of filling 
factors, N a p = v a 8 a p, and also write the fields and their con- 
jugates as vectors in the component space, a° = (a°, . . . , a°) 
and P° = (Pi, . . . ,P°). We omit the q dependence for a 
while. In this concise notation, the oscillator Hamiltonian can 
be written as 



l?_eB_ 
2m h 



[e 2 a°^N a° + b 2 P 0<i K^NK P°] . (3.8) 



The prefactor can also be written as L 2 ui c /2h, where w c = 
eB/m is the cyclotron frequency. We recall that the matrix 
K is real and symmetric, so that K * = K. We perform the 
diagonalisation in two steps. First, we define a' — y/~N a° and 
P' = V N~ 1 P°, so that the Hamiltonian becomes 



L 2 lo c 
2h 



e 2 a ,<! a' + b 2 P'^VN KNKVNP' 



The matrix between the P"s is the square of the matrix E = 
^/N K\/N ', which is real and symmetric. Therefore, it can be 
diagonalised in terms of a diagonal matrix D and an orthogo- 
nal matrix C, such that E = C~ 1 D C. The matrix D has the 
eigenvalues A Q of E on its diagonal, and C T contains the cor- 
responding eigenvectors as columns. The ability to choose C 
as an orthogonal matrix (i.e., C _1 = C T ) is provided by the 
property that the matrix E is symmetric, so that the eigenvec- 
tors can be chosen such that they form an orthonormal basis. 
Having found the diagonalisation E = C T D C, we define 



Ca' = CVNa° : P = C P' = CVN^P , (3.9) 



so that the Hamiltonian becomes 

L 2 uj, 



2h 
L 2 uj ( 
2h 



\e 2 a) a - 



E 



6 2 P t £) 2 P 



e 2 ala a + b'P' a X 2 a P a 



2 E>t \2 D 



(3.10a) 
(3.10b) 



written in matrix form and in components, respectively. For 
the derivation we have used that J2 a ®ofl<x — — 
a''a' by virtue of the orthogonality of C, ^ 7 C ai Cp 7 = 
Y^-y Ca-yC^p — 8 a p. For this transformation to be well- 
defined, it is required that vp ^ for all components f3, which 
we already assumed in order for the harmonic approximation 
to be valid. The definition is such that the commutator be- 
tween a and P is given by 



[a Q (q),P /3 (-q')] = ift<5 Q/ 3<5 q , q ', 



(3.11) 



which holds also by virtue of the orthogonality of C. 
By setting P° = — ifr-J^r and consequently P a 



-iK-sMr-, we can derive that 



= exp^-^££a a (-q)£ Q a Q (q)^ (3.12) 



is a ground state of the Hamiltonian ( 13.101 ) if we set £ Q = 
|A Q |~ . Evidently, £ Q is only well defined if the matrix E is 
nonsingular, i.e., if none of its eigenvalues is zero. Moreover, 
the eigenvalues appearing in the eigenstate are actually not the 
eigenvalues of E itself, but the square roots of the eigenvalues 
of E 2 = y/ NKNKy/N, namely ^/X 2 = \X a \. The ground 
state ( 13.121 ) can then be written in matrix form as 



exp 



2Kb 



exp 



2hb 



(3.13) 



where we used a' D^ 1 !! = a°t K~ x a° in order to write the 
ground state in terms of the original auxiliary gauge fields 
a° . Notice that, had we chosen the negative square roots 
— y/X 2 ^ for the eigenvalues of E, the ground-state wave func- 
tion ( 13.13b could not be normalised. Negative eigenvalues are 
indeed unphysical because they would lead to an instability of 
the electron liquid, the components of which phase-separate, 
as may be seen within the plasma picture of the FQHEi 2 ^ It 
is therefore important, for the structure of the Chern-Simons 
theory to be well-defined, to discard negative eigenvalues A Q . 
This is namely the case for the analysis presented in Sec. lIIICl 
where we assume a positive definite K, The case of zero 
eigenvalues is treated separately in Sec.lIVI 

Acting with the Hamiltonian ( 13.101 ) on the ground state 
( 13.12b gives its energy eigenvalues 



^Ew = ?E< 



where uj a — |A Q |o; c are the characteristic frequencies, given 
in terms of the eigenvalues X a and the cyclotron frequency uj c . 
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At this point, we define the ladder operators as 

•Aa{<ti = n (e a Q (q) + ib\ a P a (qj) , 

^ _ (3.14) 

At (<j) = ^ 47r ^ A (e3«(-q) - i6A a P a (-q)) , 

still under the assumption that the eigenvalues A Q are pos- 
itive. The commutator of the rescaled ladder operators be- 
comes [^ Q (q),^(q')] = SapSqrf, so that .4+ (q)An(q) is 
the number operator for the oscillator states in the component 
a of the diagonalised basis. The Hamiltonian can be conve- 
niently written in terms of the ladder operators as 

ffo S c = ^^^(ii(q)i«(q) + |). (3.15) 

q a 

This result also proves that the "ground state" (13.12b is indeed 
the lowest-energy state. 

Notice that the energies fnu a play the role of quasi-p article 
gaps in the Chern-Simons theory, and the ground state is well- 
defined for dct(K) ^ 0. 8 Zero -energy gaps are obtained if 
one of the eigenvalues X a = 0, i.e., when the matrix K is 
singular, det(K) = det(E) = 0. Contrary to what one may 
naively expect, this situation is not in contradiction with an 
incompressible quantum liquid, where all (collective) charge 
modes must be gapped. As we discuss in more detail in 
Sec. IIVI the zero-gap modes associated with A Q = re- 
veal ferromagnetic properties of the underlying state, 22 which 
in the presence of interactions evolve into spin-wave modes 
while keeping the charge modes gapped. 



C. Connection with trial wave functions 



In order to obtain the wave functions corresponding to the 
ground state ( 13.131 ). we may rewrite it in terms of the density 
fluctuations Sp a (q), using the constraint ( 13.31 ). Once again, 
it is more convenient to do the computation in matrix nota- 
tion. The constraint is then given by a = (h/e\q\)K (8p) for 
physical states, with (8p) = (Spi, . . . , 5p K ) the vector of the 
density fluctuations. Hence, we find 




(3.16) 



Notice that, written in terms of density fluctuations, the 
ground-state wave function is no longer confronted with the 
problem of zero-eigenvalues of E (or K) because it is the ma- 
trix K, and not its inverse K~ l , which appears here. 

As shown in Ref. \2l[ we may relate the expression ( 13.161 ) 
to the plasma picture proposed by Laughlin in his original 
publication. 1 In this picture, we regard |xosc| 2 as the Boltz- 
mann weight exp(— f3W) of the plasma Hamiltonian H, where 
one sets (3 = 2 (Ref. |26J). Then T~L can be identified as the 
Hamiltonian of particles interacting due to the Coulomb po- 
tential in two dimensions, — log | r | , which equals 2ttL 2 /\q\ 2 
in momentum space. As discussed in Appendix [A] the wave 
function that we obtain is 



^(Kj)=n n 



z k n 



Ot j a .k a 

j a <k a 



n IK 

a,p ja,kp 



Zk B ) Kaf> exp ( - ^ v a K a p ^ 

a,fi kp 



\ z ka 



<f>{»-}({ziJ), (3-17) 



where we write z — x — iy. This wave function is a product 
of the oscillator function and the wave function (f>r v » j ({Zj a }), 
which encodes the residual degrees of freedom for particles 
in the reduced field B Q , i.e., at the effective filling factors 
v* a given by Eq. (12.5b . Quite generally, one may describe the 
same system in the framework of different Chern-Simons the- 
ories, according to how much flux is absorbed in the transfor- 
mation by the matrix K a p. It is often convenient, if possi- 
ble, to choose the Chern-Simons transformation such that the 
residual wave function is factorisable into single-component 
wave functions <f> v * , 

K 

^}({^j)=n^s(^»' < ii8) 

a = l 

so that each component may be treated independently after the 
transformation. Notice, however, that this aim may be in con- 



flict with the above-mentioned condition of positive eigenval- 
ues of the charge matrix K a p, namely in the context of sym- 
metric states with ferromagnetic properties that we discuss in 
Sec. UVBl 

The simplest state of a factorisable residual wave function 
according to Eq. ( 13.181 ) consists of a product of states at an 
effective filling factor v* a = 1 for each component, each of 
which involves a Slater determinant, in the form 

j a <k a ^ k a B ' a ' 

(3.19) 

Such a state would then correspond to a Halperin wave func- 
tion that is described by an exponent matrix M a p = K a p + 
8 a p. In order to have a fermionic wave function, the elements 
K aa must naturally be even integers, and we thus have a 
Chern-Simons theory that transforms fermions into (compos- 
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ite) fermions. Alternatively, one may have chosen the bosonic 
version of the Chern-Simons theory, in which case the diag- 
onal elements of the matrix K a p = M a p would be odd. The 
same state ( 13.171 ) would then be described as a product of the 
oscillator wave function Xosc, which absorbs all the flux, and 
a bosonic wave function for zero net magnetic field £>* = 0, 
for all components, = o} ({ z j a }) = 1- 

Until now, we have discussed states that may be de- 
scribed in terms of generalised K-component Halperin wave 
functions, 22 where the residual wave function (/>{„* y {{zj a }) is 
itself a (typically simpler) Halperin wave function described 
by a "residual" exponent matrix M*g such that M a p — 
K a p + M*g (see also Appendix [A). Notice, however, that 
the Chern-Simons theory discussed above may also provide 
us with another class of factorisable trial wave functions if 
we replace the Slater determinants ( 13.191 l for the effective 
filling factors v* a = 1 by Slater determinants for p a com- 
pletely filled composite-fermion levels <f>p] ({%,}) in each 
component. The resulting wave function ( 13.17b is related to 
the re-component Halperin wave function in the same manner 
as Jain's one-component composite-fermion 2 - 3 to Laughlin's 
wave function. 1 Naturally, the proposed Slater determinants 
contain non-analytic components in the polynomial, and, in 
the same manner as for Jain's wave functions, one needs to 
project the resulting wave function to the subspace of analytic 
functions in order to satisfy the lowest-Landau-level condi- 
tion. 

Ultimately, the theory may be generalised to the case where 
the i/*'s can take any fractional value, as to allow the multi- 
component generalisation of higher-generation FQHE states. 
An example of the latter in one component is the v = 4/11 
state, which can be understood as a second generation FQHE 
stated 



IV. SINGULAR TRANSFORMATIONS 

The analysis in the previous section demonstrates that a 
Chern-Simons transformation with a nonsingular charge ma- 
trix is already interesting in itself. However, transformations 
with singular charge matrices play an important role in the 
study of states with (partial) ferromagnetic order, since these 
states are described by singular exponent matrices. 22 In this 
section, we investigate the consequences of the symmetry 
properties of the exponent matrices M and M* and the charge 
matrix K for the results of the previous section. 

A. Conditions on the ranks of the matrices 

Without performing the diagonalisation of the oscillator 
Hamiltonian, it is akeady possible to give some conditions 
on the exponent matrices and the charge matrix. Con- 
sider a state that is described by a singular exponent ma- 
trix M. As a consequence, not all filling factors are de- 
fined separately. Suppose furthermore that the electronic and 
composite-fermion filling factors are given by Y^p M a pvp = 
1 and J2b m Ib v I = respectively, with M = M* + K. We 



note that Eq. ( 12.51 l has to be satisfied simultaneously, which 
does not necessarily follow from the other conditions. 30 From 
the fact that M, M* and K are required to be nonnegative 
definite, it follows that also K and M* are singular. More 
specifically, it follows that the null spaces of M* and K may 
be of higher dimension than that of M. As a consequence, the 
dimension of the null space of the exponent matrix is either 
increased or kept invariant by the Chern-Simons transforma- 
tion. In other words, if before applying the Chern-Simons 
transformation the theory involves a certain number of inde- 
pendent combinations of filling factors, then the number of in- 
dependent combinations after the transformation is either the 
same or lower. In terms of the ranks of the matrices, which is 
equal to their size minus the dimension of the null space (i.e., 
dim ker M + rank M — k), we find that the ranks of K and 
M* must both be smaller than or equal to the rank of M. 

For the case that rankM* < rank M, which is not ruled 
out by the above discussion, some problems may arise. In this 
case, Eq. ( 12.51 ) fixes the filling factors v* a to be confined to a 
subspace of the space of all solutions of M*gVp = 1- 
For example, if M = (13) and K = ( q P, ) , we have 
M* = ( I I ), so that, based on the exponent matrices, the 
electronic and composite-fermion filling factors are given by 
v 2 ) = (1/4, 1/4) and v\-\-v\ = 1, respectively. How- 
ever, based on Eq. ( 12.51 ). the composite-fermion filling fac- 
tors are fixed at (z/f, — (1/2, 1/2). Therefore, the matrix 
M* does not appropriately describe the possible composite- 
fermion filling factors of the system. We would expect that 
this leads to problematic results, if we used the Chern-Simons 
approach to obtain a separation between high-energy and low- 
energy degrees of freedom. For this reason, we will only anal- 
yse the case that M and M* share their ranks. We stress that 
there is no problem in using a singular charge matrix K if M 
and M* are both nonsingular. 

B. The oscillator Hamiltonian 

Here, we discuss how the singularity of the matrix K af- 
fects the analysis that we used to study the harmonic oscilla- 
tor. Apart from the zero modes in the harmonic oscillator, we 
must also take into account that the number of independent 
constraints [Eq. (13. 3H is reduced, since these also involve the 
matrix K. Indeed, the number of independent constraints is 
given by the rank r of the matrix K, whereas the number of 
zero modes is k — r. Before we derive the fully general re- 
sults, we find it instructive to illustrate the procedure first with 
a simple example. 

We consider a two-component system, where we choose 
the charge matrix of the Chern-Simons transformation to be 
the singular matrix K — (22)- The eigenvalues of K are 
4 and 0, and the respective eigenvectors are (1, l)/v2 and 
(1, — 1)/V2. We can write the constraints in components as 

0= (a°(q)-^(2<5 Pl (q)+2^ 2 (q))) |0 phys ) , 

for a = 1,2. The two components a\ and of the gauge 
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field satisfy the same constraint, so that they are fixed to the 
fluctuations of the total density Spi + 8p 2 . On the other hand, 
the difference of density fluctuations Spi — 5p 2 (associated 
with the zero eigenvalue) is absent, implying that one may 
have zero-energy fluctuations that lower the particle number 
in one component while increasing that in the other compo- 
nent. Eventually such a reorganisation of the particles on the 
two components may even completely polarise the system, 
with v± = v and v 2 — 0. Inversely this means that in the 
case of a singular matrix K, we may always choose both fill- 
ing factors nonzero or even equal, i.e., N nonsingular, such as 
to render the harmonic approximation ( 13.7b valid. 

We now turn to the harmonic oscillator Hamiltonian. We 
write N — diag(z/i, 1/2), where 1/^2 are the electronic filling 
factors. For this example, we compute 

e = VnkVn = 2 ( * 

This matrix is diagonalised as C T DC, where D = 
diag(Ai, A2) is the diagonal matrix with the eigenvalues Ai = 
2(i^i + 1/2) and A2 = 0. The corresponding eigenvectors are 
proportional to (y/vi, \Jv2) and (—1/1^2,1/^1), respectively. 
The diagonalised Hamiltonian is given by Eq. ( |3.10bt , where 
a = 1 , 2 and 

fai\ _ 1 ( v\a\ + v 2 a° 2 \ 
\a 2 J ~ y/v x + v 2 \\/v\V2 + 0,2) J ' 

(TA 1 + \ 

{P2J V^TV 2 {-^Pi° + sJf 2 PZ) ■ 

We note that P 2 is not present in the Hamiltonian since the 

term P 2 X 2 P 2 vanishes due to A2 = 0. The term a\a 2 also 
vanishes, since a 2 = 0, due to the constraint a° = a 2 . In the 
end, we obtain a harmonic oscillator Hamiltonian with only 
one coordinate (at) and one momentum (Pi) component. 

The Hamiltonian restricted to this single coordinate has a 
ground state Xosc.i = exp[— (e/2ttb)a'D a], where we define 
D = diag(i(^i+^2) _1 7 0). We note that Xosc.i only involves 
01, but not a 2 . Transforming back to the coordinates (a°, a 2 ) 
and imposing the constraints a° = a 2 = (h/e\q\)(2Spi + 
2 5p 2 ), we obtain 

Xosc.i = exp(~(5pi + 8p 2 ) f ^^-(2) (Spi + 8p 2 )^j , 

where the notation (2) is to point out that it should 
be interpreted as a matrix. At this point, we observe 
that (&pi + Sp 2 y (2)(Spi + Sp 2 ) is exactly equal to 
(Spi,Sp2yK(Spi,Sp2). This means that in this example 
Eq. ( 13.161 ) is valid without change, and the other results con- 
cerning the Halperin wave functions hold as well, as we have 
already mentioned in the discussion of the general oscillator 
function ( 13.16b . We remark that the linear combination of fill- 
ing factors v\ — v 2 is not present at all in the diagonalised 
theory. 

Another important point is that we can make the connec- 
tion with ferromagnetic Laughlin states in two-component 



systems. 26 For instance, the exponent matrix M = ( f f ) de- 
fines a state for which the total filling factor is vx + V2 = 1 /3, 
but the separate filling factors are not defined, since the expo- 
nent matrix is singular. 26 Using the Chern-Simons transfor- 
mation of the example above, we may understand this state 
in terms of a composite-fermion theory with exponent matrix 
AI* = (1 !)■ This state has total composite-fermion filling 
factor v\ + v 2 = 1, and again the separate filling factors are 
undefined. We remark that although the intermediate steps in 
the procedure contain the separate filling factors v\ and v 2 , 
the results are completely independent of v\ — v 2 . 

In contrast to the ferromagnetic Laughlin state discussed 
in the preceding paragraph, we may also discuss the two- 
component state at total filling factor v = 2/5, described 
by the matrix M — (23) ar, d me reduced exponent matrix 
AI* = (oi). Although the Chern-Simons transformation 
is described by a singular charge matrix K and does there- 
fore not impose a constraint on the relative particle distri- 
bution on the two components, the constraint is imposed by 
AI*, v\ = v 2 = 1. The state thus described is then a spin- 
unpolarised state, as one could have also expected from the 
original exponent matrix M, 

The reasoning given for the example above can be read- 
ily generalised to any situation in which K is singular. Here, 
we assume that the rank r of the matrix K is smaller than 
the number of components k. As argued in Appendix iBl the 
ground state can be decomposed as a product of the usual 
ground state ( 13.131 ) restricted to the r independent compo- 
nents, Xosc.r and the degenerate part \ t see Eqs. ( IBU and 
dB2M . Moreover, Eq. ( 13.16b remains valid even in the singular 
case, despite the fact that the original derivation involves the 
inverse of K. Hence, the Halperin connection in Sec. IIII Cl is 
valid in the singular case without modification. 

The equivalence of the decomposition ( IBU for the n- 
component oscillator wave function may be interpreted in a 
straightforward physical manner. Indeed, the decomposition 
indicates that, in the case of a charge matrix K of rank r, 
the "reduced" r-component wave function corresponds to an 
r-component Halperin wave function with gapped oscillator 
frequencies u a . The other factor \ m Eq. (IB II) corresponds 
to the k — r zero eigenvalues of the matrix K with an as- 
sociated space spanned by the oscillator components a a , with 
a = r + 1, . . . , K. The ground-state manifold comprises there- 
fore any possible combination of these components a a , and a 
particular choice spontaneously breaks the residual ground- 
state symmetry, which may be related to the ferromagnetic 
properties of the Halperin state, and \ may then be interpreted 
as the ferromagnetic part of the wave function. 

In order to see this particular point, consider the r con- 
straints to fix the filling factors of the first r — 1 components. 
The last constraint then imposes simply the sum of the fill- 
ings of all other components a = r, . . . , k. This is naturally 
a simplified assumption, because the r constraints do not in 
general fix particular components, but the dependencies may 
be more complicated. 31 One is then free to distribute the in- 
volved particles over these components in a quantum mechan- 
ical manner. All different distributions define the ground-state 
manifold. Schematically, this may be formalised with the help 
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of a wave function 

X = u r \a — r) + u r+ i \a = r + 1) + . . . + u K \a = k) , 

where the complex amplitudes u a are subject to a normali- 
sation condition, which plays the role of the last constraint. 
These complex amplitudes may be viewed as the components 
of a CP K ~ r fields The ground-state manifold may then be 
described by spatially constant CP K ~ r fields with a global 
SU(k — r + 1) symmetry, which is precisely the symmetry 
group that describes the ferromagnetic properties of the oscil- 
lator wave function. In summary, this argument shows that, 
in the case of a Chern-Simons transformation with a matrix 
K of rank r, one may decompose an arbitrary oscillator wave 
function into a product of a reduced r-component Halperin 
wave function and a S\J(n — r + l)-symmetric ferromagnetic 
part. Naturally, this symmetry may be further reduced if the 
components of the Chern-Simons field fix further filling fac- 
tors. 

We finally mention that the spontaneous breaking of the 
SU(k — r + 1) symmetry yields Goldstone modes, which are 
physical (pseudo)spin waves. On the level of the Gaussian 
model, these Goldstone modes are dispersionless and remain 
at zero energy. This is no longer the case if one takes into 
account interactions between the particles associated with the 
different components. One may indeed treat rather easily a 
density-density interaction within the present model. This in- 
teraction may be translated, via the constraints (13.3b , into an 
interaction between the oscillator fields, which one can then 
diagonalise within the Gaussian model. Notice that these 
fields are coupled to the U a [see Eq. d3.5H , which describe 
the low-energy electronic degrees of freedom. A discussion 
of collective Goldstone-type modes is therefore more involved 
and requires a decoupling of the oscillator and the electronic 
degrees of freedom. However, the Chern-Simons analysis 
within the Gaussian model yields valuable insight into the 
ferromagnetic properties of the states, which are governed 
by symmetry, as well as into the number of their Goldstone 
modes. 



saliently, one needs to discard charge matrices with nega- 
tive eigenvalues because the associated Chern-Simons theo- 
ries yield oscillator ground-state wave functions that are not 
normalised. This is in line with physical insight obtained from 
a multi-component version of Laughlin's plasma picture^ ac- 
cording to which charge matrices with negative eigenvalues 
yield inhomogeneous ground states where the components 
phase-separated 

Whereas singular charge matrices, with zero eigenvalues, 
had originally been discussed by Lopez and Fradkin 10 only for 
the SU(2)-symmetric case, we have argued here that the as- 
sociated Chern-Simons theories reflect underlying ferromag- 
netic states in a more general setting. Indeed, we have shown 
that the density fluctuations of the n components are then de- 
termined by only r < n constraints, such that n — r partic- 
ular combinations of the component densities may be chosen 
freely in the ground-state manifold, which is thus described 
by the SU(k — r + 1) group. This symmetry is spontaneously 
broken by a particular ferromagnetic state, which can be de- 
scribed by k — v different Goldstone modes that may be viewed 
as generalised spin waves. Our results encompass the particu- 
lar SU(2) case of two-component Chern-Simons theories dis- 
cussed in the literature^— 

We emphasise moreover that the analysis of the micro- 
scopic multi-component Chern-Simons theory within the 
Gaussian approximation heuristically yields trial wave func- 
tions for multi-component quantum Hall systems that may 
be further studied numerically. As an example, we have dis- 
cussed generalised /{-component Halperin wave functions that 
play a similarly central role as Laughlin's wave functions do 
in one-component quantum Hall systems. Beyond these gen- 
eralised Halperin wave functions, we have briefly discussed a 
second class of states, where the residual wave function that 
is not encoded in the Chern-Simons oscillator part % osc is 
a product of Slater determinants of p a completely filled (a- 
component) composite-fermion levels. This construction is 
reminiscent of Jain's generalisation of one-component Laugh- 
lin wave functions to filling factors v = p/ (2sp + 



V. CONCLUSIONS 

In conclusion, we have studied a microscopic Chern- 
Simons approach to general multi-component quantum Hall 
systems (with k components). Beyond the mean-field ap- 
proximation, which yields a renormalisation of the magnetic 
field that depends on the average particle densities for each 
component, their fluctuations are taken into account within a 
Gaussian model of auxiliary gauge fields. These gauge fields, 
introduced by Shankar and Murthy in the framework of the 
Hamiltonian theory of the FQHEr^— are indeed connected 
via constraints to the component density fluctuations. 

The analysis of the Gaussian model — although it may 
be viewed as a first step in the discussion of a more com- 
plete Hamiltonian theory for multi-component quantum Hall 
systems — already yields valuable insight into the structure 
and the correctness of the Chern-Simons theory, which is 
characterised by a symmetric n x k charge matrix K. 9 Most 
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Appendix A: Multi-component plasma analogy 

The single-component plasma analogy proposed by 
Laughlin 1 is readily generalised to the multi-component case. 
Here, we use the ground state 

/ 1^^ , n 2tt£ 2 . , \ 
Xosc = expl -- ^^dp^-cti—^K^Sp^q) I, 

(Al) 
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which is Eq. (I3.16l l written out in components. Recalling that fluctuations Sp a (r) — J^j <K r ~ r j a ) ~ n a- Then, we can 
27ri 2 /|q| 2 is the Fourier transform of — log|r|, we perform an rewrite Xosc as 
inverse Fourier transformation and we substitute the density 



exp 



~ K aP [ d 2 rd 2 r' ^(r-rJ-J log|r - r'| ( £ 6(r> - r fc , ) 



By evaluating the integrals, one finds 



X osc = const • Yl Yl \ rja -r k? \ K ^/ 2 exp(-^^2n a K a ^\r kfi \ 2 ] 

a, 13 j a ,k/j ^ a, j3 kp 

jc^kf, 



Using that v a — 2Trl 2 B n a , and changing to complex notation, with z = x — iy, 33 we can explicitly write this expression as 

2 - 



Xo 



const • Y\ Yl 



% Q - z k c 



a ja,k a 

j a <k a 



n n 

a, 13 j a ,kp 
a</3 



Z ja Z k, 



exp( - v a K a i3 

a, /3 ka 



\ z k a 



4/ 



B 



The Jastrow-like products in this expression only contain distances between the particles, i.e., only the moduli \zj 



Zk B 



Phase factors of the form [{zj a — z kfs )/\zj a — z k „ 



exp[iK a p arg(z Ja - z kfj )] = exp[-iK a/3 9(r ja - r kfj )} are obtained 



from substitution of the full density p a (r) — J^j ^( r — r j a ) mt0 the Chern-Simons transformation d2.1| ). Applying this 
transformation to Xosc, we obtain the product of the latter with the phase factors, 

\ r h„ B / 



a 3c,k 
j a <k 



a,/3 j a ,kp 
a<!3 



where j denotes the composite-particle wave function for filling factors v* a , which will be investigated in the following. The 
magnetic lengths appearing in (/>{„» } ({ z j a }) are the reduced magnetic lengths Ib* given by Eq. (12.51 . 

As an example, we consider the situation in which v* a can be determined by an exponent matrix M* , 9 ' 22,26 such that (f>{„* j is 
the Halperin wave function 



*>({*i.})=n n fe^o^n n^-^^f-EE 



a j a ,k a 
j a <k a 



Q</3 



4/ 2 



Combining Eqs. (IA21 > and dA3l ), we obtain the full electronic wave function, 



(A3) 



v>({zij)=n n ( 



a ja-k a 
j a <k a 



Ot,P ja,kp 

a</3 



\ z k 



4/ 



(A4) 



which is the Halperin wave function for the exponent matrix 
Ma/3 = Ma/3 + K a ^M Here, we have expressed the effective 
magnetic lengths in the exponential of Eq. ( IA3b in terms of 
the original one, as 1/1%, = (1 — Y^,a K a pvp) /l 2 B , by virtue 
ofEq. (O. 



Appendix B: Ground state in the singular case 

The reasoning given for the two-component example in 
Sec. lIV Bl can be extended to any number of components. Sup- 
pose that the charge matrix K (being a k x k symmetric non- 
negative definite matrix) is of rank r, which means that it has 
r independent rows or columns. In particular, there are k — r 
rows or columns that can be written as a linear combination 
of the other r independent rows or columns. This also means 
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that the dimension of the null space, or equivalently, the mul- 
tiplicity of zero eigenvalues is equal to k — r. 

Since the constraints ( 13.31 ) are expressed as a linear relation 
involving the matrix K, there are only r independent con- 
straints. Hence, the vector a° = (a°, . . . , a° ) lives only in 
an r-dimensional subspace; n — r of its components can be 
written as a linear combination of the other r. 

Now we analyse the Hamiltonian (13.7b . Since we have as- 
sumed that the matrix of densities N is nonsingular (i.e., all 
filling factors are nonzero, as required for the harmonic ap- 
proximation to be valid), the rank of E = y/~N K^/N is 
equal to the rank of K. This means that E has r positive 
eigenvalues and k — r zero eigenvalues, just as the matrix 
K. We diagonalise E as usual in terms of a diagonal ma- 
trix D and an orthogonal matrix C such that E = C T D C. 
Note that the order of the eigenvalues on the diagonal of 
D (and simultaneously the order of the rows of C) may be 
chosen at will, so that we may choose for simplicity D = 
diag(Ai , . . . , X r , 0, . . . , 0), where Ai , . . . , A r are the positive 
eigenvalues of E. In the diagonalised Hamiltonian ( 13.10b . the 
components P r +i, • • • , P/% are absent since they are multi- 
plied with the zero eigenvalues of D. We still have k com- 
ponents of a in the Hamiltonian, but we should remember that 
only r of them are independent. 

The diagonalised Hamiltonian contains r nonzero eigenval- 
ues, which depend on the filling factors v a . However, some 
variations of the filling factors will leave the eigenvalues, and 
hence the diagonalised Hamiltonian, invariant, namely those 
satisfying the equation 

= (VXp) ■ 8v = V -^-8v a for all p. 

a 

In other words, the desired variations are the vectors in the null 
space of the gradient matrix (VA) of the eigenvalues, which 
is defined as the matrix of derivatives of A, with respect to 
v, (VA) ) a a = dXp/dv a . Since k — r of the eigenvalues A^ 
are zero, the rank of the gradient matrix is at most r, and this 
consequently means that we can find at least k— r independent 
variations of the filling factors which leave the eigenvalues 
invariant. 

In the example discussed in Sec. IIV Bl we observed that 
v\ — v 2 does not appear in the diagonalised Hamiltonian. In 
order to demonstrate the procedure sketched in the preceding 
paragraph, we compute the gradients of the eigenvalues and 
2{v\ + v-x). Obviously, in this example the gradient matrix 
is (VA) = (22) ar, d i ts nu H s P ace is spanned by the single 
vector (1,-1). Since this vector is independent of the fill- 
ing factors v\ and v 2 , we can state that all eigenvalues, and 



hence the diagonalised Hamiltonian, are invariant under the 
transformation {v\ — > Vy + 8v, v 2 — > v 2 — 8v}. This means 
that the linear combination v\ — v 2 is completely absent from 
the Hamiltonian, as argued earlier. We remark that the varia- 
tion that leaves the Hamiltonian invariant need not always be 
constant in the filling factors v a . However, for physically rel- 
evant systems, the variations are constant, describing particle 
exchange among different components. 

We now return to the diagonalised Hamiltonian, and try to 
find the lowest-energy states, in the same way as we have done 
for the example in Sec. IIV B| For the moment, we do not im- 
pose the constraints, thus regarding all components a a as in- 
dependent. Only the first r components a„ (a = 1 , . . . , r) 
have a corresponding momentum operator in the Hamil- 
tonian, while the other k — r components do not. This 
means that the resulting states are degenerate in the coordi- 
nates a r +i, ■ ■ ■ ,a K . Thus, we may write the lowest-energy 
states as 

Xo SC (ai, ...,a K ) = Xosc,r(ai, . . . ,a r ) x(a r+ i, . . . ,a K ), 

(Bl) 

where x is me degenerate part of the wave function (further 
discussed in Sec. IIV Bl . and 

Xosc,r(ai, • ■ • ,a r ) = expf ^alX^aA (B2) 

^ a — l ' 

is the nondegenerate part. Notice that the components as- 
sociated with the zero eigenvalues of the matrix E do not 
contribute. By the observation that the pseudoinverse 28 
of D — diag(Ai, . . . , A r , 0, . . . , 0) is equal to D = 
diag(Af , . . . , A" 1 , 0, . . . , 0), we may also rewrite Xosc.r as 

Xo SC ,r = exp(-^atBa) = exp(-^a°t£a°) , 

where we used that a) Da — a°^Ka°. This result is nothing 
else than Eq. ( 13.131 ) with the inverses of D and K replaced 
by the pseudoinverses. Substituting the density fluctuations 
Sp a for the gauge fields a° using the constraint (13.31 ) yields 
exactly ( 13.161 ) by virtue of the property of K that K K K — 
K. This result is exactly equal to the steps we followed before, 
but only with K replaced by K. Therefore, the ground 
state ( 13.16b found for the case of strictly positive eigenvalues 
is also valid if there are zero eigenvalues. All subsequent steps 
concerning the connection to the trial wave functions remain 
valid as well. 
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